function val=f(x)

val = zeros(size(x));

val(x > 0.008856) = x(x > 0.008856).^(1/3);
val(x <= 0.008856) = 7.787.*x(x <= 0.008856) + 16/116;